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Abstract 



Let the adiabatic invariant of action variable in slow-fast Hamiltonian system with 
two degrees of freedom have two limiting values along the trajectories as time tends 
to infinity. The difference of two limits is exponentially small in analytic systems. 
An iso-energetic reduction and canonical transformations are applied to transform the 
slow-fast systems to form of systems depending on slowly varying parameters in a 
complexified phase space. On the basis of this method an estimate for the accuracy of 
conservation of adiabatic invariant is given for such systems. 



1 Introduction 

Consider a Hamiltonian system with two degrees of freedom. The Hamiltonian E depends 
slowly on coordinate e~^x and fast on coordinate q: 

E^E{p,q,y,x) 

where g, e~^x are coordinates, and p, y are their conjugate momenta. 

The variation of the variables p, g, y, x is described by the differential equations 



oq op 



dE . dE 



dE 
dx ' 



dE 



y 



X 



dy ■ 
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The variables p, q are called fast, and x are called slow variables. This system is called a 
slow-fast Hamiltonian system. The system with one degree of freedom in which y = const, 
X = const is called unperturbed or fast system. Assume that when y,x = const, the phase 
portrait of the unperturbed system contains a domain filled by closed trajectories [T]: 



i 




Figure 1 

Here the frequency of motion is non-zero. So we can introduce action-angle variables 

I = I{p,<l:y:X), if = ip{p,q,y,x) mod 27r 

The action / represents the area surrounded by each trajectory. It is equal to this area 
divided by 27r. Denote HQ{I,y,x) the Hamiltonian E expressed in terms of J, y, x. The 
approximation in which / = const and dynamics of y, x is described by Hamiltonian system 
with Hamiltonian HQ{I,y,x) is called an adiabatic approximation. In this approximation, 

E{p,q,y,x) = Ho{I,y,x), 

dHo{I,y,x) . dHo{I,y,x) 

y = -e 7, , x = e . 

ox ay 

We will assume that in the phase portrait of this system there is a domain in which along 
trajectories — )■ oo, ?/ — )■ const as time t ±00. Thus trajectories have a form shown in 
one of figures below: 
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(1) (2) (3) 



Figure 2 

We consider the case when the adiabatic invariant of the action / along a trajectory of the 
original system possesses limiting values /-t as t — )■ ±oo. Their difference A/ = /+ — /_ is 
called an accuracy of persistence of the adiabatic invariant in slow-fast systems. 

If this system is analytic, the value of A/ is exponentially small pQ: 

A/ = 0(e"?), 7 = const > 0. 

The goal of this paper is to find an estimate of constant 7. After constructing some assump- 
tions, we will then formulate and prove the conclusion. 



2 Reduction to a standard form 

For the original Hamiltonian E{p,q,y,x) with y,x = const, q is the coordinate and p is 
conjugate momentum. In the unperturbed system, denote T the period of the trajectory and 
u the frequency. Let S{I, q, y, x) denote the generating function of canonical transformation 
(p, g) {I, if). We have the relations [3]: 

dS{I,q,y,x) dS{I,q,y,x) 

^ = — m — ' ^ = — dq — ■ 

For a fixed trajectory E{p,q,y,x) = h = Hq{I ,y,x), where /i is a constant, we can express 
p as 

p = P{h,q,y,x) 
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dE 

if — — 7^ 0. So the generating function has the following form R]: 
op 



S{I, q, y,x) = j P{h{I, y, x), qi, y, x) dqi, 

<?0 

where the constant go is the initial value of qi. 
Lemma 2.1. (see, e.g. fj]) 



dS ff/dE\ dE\ 



dx J \ \ dx / dx J 



dti, 







dS }f/dE\ dE 



9y J \\ 9y / dy 



dti 







where 

'dE\ 1 r dE , /dE\ 1 r dE 



dx I T J dx \dy / T J dy 

Now consider another canonical transformation (p, q, y, x) t— )■ (/, ip, y, x) with generating 
function e~^yx + S{I,q,y,x). The old conjugate variables are (p, g) and {y,e~^x), and new 
pairs after transformation are (/, and {y,e~^x). 

Lemma 2.2. (see, e.g. ^l]) 

y = y + 0{e), x = x + 0{e), I = I + 0{e), ip = ^ + 0{e) 
and new Hamiltonian is 

H = Ho{I,y,x) + eHi{I,ip,y,x,e). 

We will consider the case that there exist limiting values of action variables / and I as time 
tends to infinity. The conditions for this and the proof of existence of limits will be added 
later. Then we will see that the limiting values of / and / are equal. 

So far we have got a new Hamiltonian 

Hil,(p,y,x,e) = Ho{I,y,x) + eHi{I,(p,y,x,e). 

For simplicity, omitting the bar symbols and dependence on e, we obtain a Hamiltonian 
system in standard form: 

Hil,^,y,x) = Ho{I,y,x) + eHi{I,ip,y,x) 
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which has motion 



• dHi dHo dHi 

^ ^ \ dx ^ dx ) ' ^ \ dy dy J 



3 Statement of the problem 

It is assumed, that in adiabatic approximation, \x\ — )■ oo, ?/ — > const as t — )■ ±00 (see 
Figure 2). We will assume also, that when x, y are changing in accordance with adiabatic 
approximation, and p, q are fixed: 

dE{p,q,x,y) 

> 0, as t — )■ ±00. 

ox 

For definiteness we will consider the case when the trajectories in adiabatic approximation 
have the form shown in Figure 2(1). Result will be valid for the other two cases of Figure 2. 

The function 

dHo{I,y,x) 
(^o[I,y,x) = — 

is the frequency of unperturbed motion. 

Let us assume that the following conditions are fulfilled. 

Assumption 1°. 

The functions Hq, Hi can be analytically extended into a complex domain D = DjxD^pXDxy, 
where Dj is a neighbourhood of a given real point 




* Re/ 



Figure 3 
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is a strip of a fixed width about the real axis, 

Ttn f 

I /' 

> Re 5? 



Figure 4 




and Dxy is some domain in the complex plane C^. The function uo{I,y,x) does not vanish 
in D and |a;o| > const. Function Ho{I,y,x) satisfies 



dHo 



dl 



< const, 



dy 



< const, 



dx 



, dHoY fdHo\ 
< const, I — — I + I — — I > const. 



dy 



dx 



Now let us consider approximate Hamiltonian Ho{I,y,x) which has motion 



y = -e- 



dHo 
dx 



X = e- 



dHo 
dy 



For a fixed /q, we have HQ{lQ,y,x) = h = const. Assume that curves Hq = const are not 
closed: 




h = hf)= const 



Figure 5 



Introduce the slow time 

r = et. 

The differential equations of motion are 

dy dHo dx dHo 
dr dx ' dr dy 
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Take h = ho, /iq is a point in some interval Dh with centre h^,, where h^, is a given real point: 



Figure 6 



We can find the solution for describing the motion for approximate Hamiltonian with 
Ho{Io,y,x) = ho: 

y = YiT, /o, ^o) 
X = X{T,Io,ho) 

From now on we omit the dependence on ho. 



lo e Di, ho e Dh 



Assumption 2°. 

The solutions Y{t,Io), X(r, Jq) can be analytically extended into a strip 

Dr = {t: |Im t\ < a + 6\Re r|}. 



imx 




Rer 



Figure 7 



We suppose that for any initial point, the solution tends to infinity as time tends to infinity, 
i.e. 

lim \X{r, Jo) I = oo. 



Ho{I,y,x) satisfies 



ago(/,y(r,/o),X(r,Jo)) _^ 

Rer-s-itoo dx 
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Hi (/, ip, y, x) satisfies 

|//i(7,^,y(T,7o),^(r,/o))| < 



1 + 



Ju;o{Io,Y{n,Io),X{n,Io))dn 





2+1/ ■ 



ujQ{I,y,x) satisfies 

Im uJo{Io, Y{t, /q), X{t, Jo)) =^ 0, as Re r ^ ±oo, Im Iq 0. 
Here a, 5, c, v are positive constants. 

Lemma 3.1. For I E Dj = Dj — 6j, cp e = D^p — 5^, and 
T eDr ^ {t: \lm t\ < a - 5r + 5|Re t\}, 



dHi{I,cp,Y{Tjo):X{rJo)) 



dl 



dHi{I,ip,Y{T,Io),X{T,Io)) 



< 



const 



1 + 



Juo{Io,Y{Ti,Io),X{ri,Io))dn 



2+u 1 



dip 



dHi{I,ip,Y{T,Io),X{T,Io)) 



< 



const 



1 + 



dr 



where Sj, 5^, 5^ are positive constants. 



< 



JuJo{Io:Y{riJo):X{nJo))dTi 



const 



2+u ' 



1 + 



/a;o(/o,n^i,/o),^(Ti,/o))dTi 



2+i/ ' 



Proof. Im c^o(-^o, Y{t, Jq), X(r, Iq)) ^0, as Re r ±oo, Im Jq -> 

implies that Vp > 0, 3 S > 0, F > 0, such that if |Im /o| < S, |Re r| > T, then 



|Im a;o| < p. 



Cl 



Since we know that |a;o| > Ci = const, assume |Im uio\ < — for |Im /o| < E, |Re r| > F, and 
therefore, 

9 

|Re u;o\ > —Cl. 



Let T — X + i/i. The integration 
path. 



/a;o(/o,n^i,/o),^(Ti,/o))dTi 





does not depend on the 
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For |Re r| >max|r, 2^— ci — ^—ci + /9 + a + 6^r|, where a, 6 are positive constants 
satisfying 

r r 
J |Re uJo\ dA < a ■ r, J |Im uJo\ dA < 6 ■ F, 





-1/ 9 



Im T 



** ► 



r Rer 



Figure 8 



we have 



Uo dri 





Rer 



Uq dri + J ujq dri 

Rer 





Rci 





Rc 



Rei 

Wo dA + y Wo d(i/i) 

Rer 

r Re r 



y Re Wo dA + y i Im Wo dA + y Re Wo d(i/i) + y i Im wo d(i/i) 









Rer 



Rei 

r 



Re -J 



r Rer r Rer 

Re Wo dA + y Re Wo dA + J ilm uq dX + J i Im wo dA 
r r 



+ J Re Wo d(i/i) + J i Im wo d(i/i) 

Re r Re r 
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> 



> 



> 



> 



as 
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Rci 



Rer 

r r 
+ J \I{je uJo\ dX + J\lmuJo\dX 



J Rea;odA+ J ReuJod{ifj,) — J |Ima;o|dA+ J |Im cuq] d(i/i) 



Rer 



Rer 
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9 



^^ci • (Re T - r) + — ci • i Im r - p ■ {Re r - F) + p ■ ilm r + aV + bV 

9 9 
— Ci ■ r C| ■ 1 

10 10 



p-T-p-r + aF + bT 

(l^'^l^' " ^'^^) " ('"'^1 + pr + ar + bv) 

{^c,-p)\r\-[^c, + p + a + by 
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ci-p \r 



ci-pj|r| >2(^— ci + p + a + 6jr. 



We can find a positive constant C2 < —Ci — p, such that 



r 

y" a;o(/o,n^i,/o),^(Ti,/o))dri 



> ( J^Ci - p ) |T| > C2|T|. 



For |Re r| < max |r, 2^^ci — (^"^^ + p + a + 6)r[>, it is evident that we can find 
another positive constant C3 such that 



r 



> C3\t\ 



So 



\H,{I,<p,Y{rJo):X{Tjo))\ < 



1 + 



Jujo{io,Y{n,lo),x{n,lo))dn 



•2+v 



< 



< 



1 + (const • |t|)2+^ 
const 



1 + |t|2+^' 
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From Cauchy estimate we can easily get that, in Dj and D<^, 



ai/i(/,y.,r(r,/o),X(r,Jo)) 



dl 



9ifi(/,y.,F(r,Jo),X(r,/o)) 



dip 



< 



< 



const 
1 + |r|2+^' 

const 



Now let us prove 



9ifi(J,<^,r(r,/o),X(r,/o)) 



dT 



const 
< , ,^ , ■ From 

1 + r 2+^^ 



= {r: |Im r| < a + 5|Re r|}, 



for simplicity, assume ~ ^^"^ then 



cr 



= {r: |Im r| < o- - — + (5|Re r|}. 



In the following figure, mark the boundary of D^- with solid lines, and that of Dj^ with dashed 
lines: 



Im T 




centre t. radius 



Figure 9 



For r G -D,-, the estimate of \HA on the circle of radius — around r is 

I II 

const 



^1 < 



1 + min |r + ^e^^P'^'' 



. If |r| > cr, then 



l^il < 



const 



1 + min 



?+-e" 
10 



< 



const 



9 

10' 



-r 



2+u ■ 
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So, because of the Cauchy estimate 



dHi 



dr 



< 



const 



a const 

< 



|2+. / 1 + \t 



2+1^ ■ 



If |r| < a, then because of the Cauchy estimate, 



dHi 



dr 



I a const 
< const — < 



10 l + |f 



2+1/ ■ 



Therefore, 



ai/i(/,(^,F(r,/o),X(r,/o)) 



dr 



< 



const 
1 + |r|2+^' 



We can find a constant C4 such that 



Then 



T 

J a;o(/o,y(Ti,/o),X(Ti,/o))dTi 



< C4 b" 



const 



1 + |t|2+^ 



< 



< 



const 



1 + ( C4 ^ 



/a;o(/o,i^(Ti,/o),X(ri,/o))dri 



2 + !/ 



const 



/a;o(/o,V'(ri,/o),X(ri,/o))dri 




2+u ■ 



Therefore, for / e Dj, (p e -D^, r e -D,-, 
9i/i(7,^,y(r,/o),^(r,/o)) 



dl 



9//i(7,(^,y(T,7o),^(r,/o)) 



< 



const 



2+i^ 



d(p 



9i/i(7,^,y(T,7o),^(r,/o)) 



< 



9t 



< 



1+ Jujo{io:Y{nJo):X{nJo))dn 



const 

1+ /a;o(/o,>^(Ti,/o),X(Ti,/o))dri 



const 

2+i7 



1 + 



/a;o(/o,>^(ri,/o),X(ri,/o))dri 
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Assumption 3°. 

The level lines 

Im j coo{I, r(ri, /), X(ri, /)) dn = B = const, < |5| < 7, I E Dj 


lie in the domain D^- and have a positive distance from the boundary of D^. 

Now consider the exact solution I{t), ip{t), y{t), x{t) of the Hamiltonian system with Hamil- 
tonian 

Hil,(p,y,x) = Ho{I,y,x) + eHi{I,(p,y,x) 

with real initial conditions /(O), ip{0), y{0), x(0) in D = Dj x x D^y at 

t = r/e = 0, and H{I{0),(f{0),y{0),x{0)) = Hq. The adiabatic invariant of the action J(t) 

satisfies 

dJ ^ ^ dHi{I,ip,y,x) 
dt dip 

and thus 



Theorem 3.1. There exist limiting values 

J±= hm I{t) 

and their difference AI = 1^ — 1^ satisfies the estimate 

AI = 0(e"?), 7 = const > 
with the constant 7 introduced in Assumption 3°. 

Remark. The method of continuous averaging [8j gives the same estimate for 7, as the 
above-stated theorem. 



4 Example 

Paper [7] gives us several examples of slow-fast Hamiltonian systems, in one of which the 
Hamiltonian is as follows: 

H{I, ip,y,x)=uI + ^ + Vo{x) + sg{ip)Vi{x). 

13 



Here x are slow variables and J, (p are fast ones. We explicitly treat the case of Vo(x) = e~^. 
Concerning Vi, we need it to be real analytic and tend to zero as x ^ +00. For simplicity 
we choose Vi = e~^. The unperturbed Hamiltonian system is 

2 2 
i/o(/, x) = + |- + V^{x) =coI+^ + e-\ 



r 



Let Ho{I,y,x) = ujI + — + e~ 



- ho = const be the total energy of the unperturbed system. 
We can draw the phase portrait of the system to describe the motion: 




Figure 10 



Take f as the slow time of motion. From the Hamiltonian Hq{I, y, x) = ujI + + ^ ^ = ho, 



we know that 



So 



dx dHn 



df dy 



1 / dx 



Thus we can express x and y via fj): 

x = x{l,fi), y = y{^,fj) 



with initial data x^ = x{0,r]) and y^ = y'2(r] — e ^°). 

After solving differential equations, we can obtain the solutions as (see [7]^ 



77) = log Q ( cosh VP [i - i^) ) , y{l r/) = tanh ^rp - 1°) • 



14 



For simplicity we take the initial value of ^ as ^° = 0. So x° = x(0, rj) — log - , y° = 0. 

Now after the canonical transformation of variables from {x,y) to (C)^)) the Hamiltonian 
Hq becomes Kq{^, fj) = uI + rj and the new Hamiltonian is 



K{i, fj,I,^) ^ujI + f} + efii, f})g{ip) 



with 



f{tfl) = V,{x{lf,))=e 



- log (^4 (cosh I) _ 



(cosh ^/r]J2^ 



The differential equations for describing the motion are 



t) — —e 



dfj 



(p = ui. 



dip 



For g{ip), we suppose that g{ip) is analytic in the strip |lm<^| < p with constant p > 0, and 
bounded: 

\g{<p)\ < 1, for \lm(p\ < p. 



For f{^,fj), we can easily prove that 



\filv)\< 



const 



and 



const 




1 + |^|2+. 



Take as a neighbourhood of a given real point and is a strip |Im</7| < p. Now let 
us determine the domain of ^ by finding singularities of /(^, fj): 



filfj) 



V 



( cosh \/r]/2 ^ 
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Points of singularities should satisfy 



+ 1=0 



= -1 



^cos Im v^2^ f + i sin Im ^/2f] = 



So cos Im = -Q-^^'i < 0, 

sinlm \/2f)f = 0. 
Therefore, Im ^2f^f = (2A; + l)7r, /c e Z, 
and Rev^l = 0. 

Thus 

^2^1= (2A; + l)7ri, A; e Z. 

The domain of fj denoted as Dfi is a complex ball around a real value f)^ . Any value fp e Df^ 
can be taken to determine the singularities. Therefore, we obtain the singularities as 



TT 



From the above conclusion, if we take ^ in the strip 



i, kez. 



D,= {i: \\mi\ < 



IT 



'2fj 



<5|Ree|}, 



where 5 is a positive constant, f{^,ri) is analytic for any rj G -D,^. 
The following figure shows the boundary of domain D^: 
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We can find the level lines 



Im J wdfi = wimf = B = const, < |-B| < 7 


lying in the domain D^, which have a positive distance from the boundary of this domain. 
So 7 has been found as any constant such that 

7 < — ; , for such definition of fin. 

According to Theorem 1, the accuracy of conservation of adiabatic invariant AI has the 
estimate: 

AJ = 0(e-?). 

This estimate coincides with that in [7], where it was obtained by a different method. 
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5 Proof of the Theorem 



The original system is described by Hamiltonian 

H{I, (f, y, x) = Ho{I, y, x) + eHi{I, cp, y, x) 

with motion in time t: 

■ dHi dHo dHi 

fdHo dHA fdHo dHi 



dx dx J ^ \ dy dy 

Air /T \ dHQ(I,y,x) , . , . , . 

As the frequency oJo[l ,y,x) — — does not vanish m our domam D, we can express 

I from H via </?, y, x and a constant of energy. This is the procedure of isoenergetic reduction 
on the energy level H{I, (f, y, x) — ho, ho — const: 

-I{y, X, (p, ho) = Fo{y, x, ho) + eFi{y, x, cp, ho). 

Here (— /) is the new Hamiltonian, is the new time, and {y, x) is in a small neighbourhood 
of Dxy, which can be written as D^y + Sxy, where 

y = y 

X = X{T,Io,ho) 

and Sxy is some positive constant. 

We do not indicate dependence of Y, X on e. The differential equations of motion are 

dy fdFo dFA dx f dFo dFi 

= Ms ' 3~=^ Ms 



Dxy^\(y,x): I rei^., heDj}, 



dip \ dx dx J dip \ dy dy 

Since we know functions Ho and Hi are analytic in D, from the implicit function theorem, 
functions Fq and Fi are also analytic in x {D^y + S^y). 

Now consider the approximate system with Hamiltonian 

Fo{y,x,ho) = -lo- 

Introduce slow time ^ = £</?, and the motion is 

dy _ _dFo dx _ dFo 
d^ dx ' d^ dy 

We can find the solution for the above approximate system: 

X = A /o, ho) 
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In what follows, we do not indicate the dependence on constant Hq. The relation between 
solutions Y{t,Io), X(r, Jq) and F(^,/o), X(^,/o) follows from the formula 

T 

e= / uo{Io,Y{n,Io),X{TuIo))dTi. 



Solution Jq), Jq) can be analytically extended into some domain with respect 
to ^. The domain will be introduced in later lemma. As we know, 

Im wo(/o) ^(t, Iq), X{t, Jo)) =^ 0, as Re r — )• ±00, Im Jq — )■ 0. 

It implies that Uq uniformly tends to a real function depending on Re/o as time tends to 
infinity and Jq tends to real. Let us denote the hmits as co±{ReIo): 

uo{Io,Y{ti,Io),X{ti,Io)) ^ u±{ReIo), as Re r ±cx), Im/o 0. 

Lemma 5.1. For Jg G Dj close to the imaginary axis, the image of the domain under 
the map r i-> ^ given by the formula 



T 

J uo{io,Y{n,lo),x{n,lo))dn 



is a strip around the real axis that contains a domain 



|Im ^1 < 0-1 + 5i|Re ^1, |Re^|>C5r 
|Im < ^2, |Re el < CsF 
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where ai, 02, ^i, C5 are positive constants and a2 = cri + Si- c^T for continuity. Straight lines 
Im^ = ±7 lie in and have a positive distance from the boundary of (i.e. 7 < a2). 

Lemma 5.2. The function Fi{y,x,ip) satisfies the estimate: 

const 



\F,{Y{^,Io),X{^,Io),f)\ < 



1 + 



2+1/ ■ 



Proof. First of all let us list two forms of Hamiltonian: 



Ho{I,y,x) +eHi{I,(f,y,x) = ho, 
-I = Fo{y,x)+eFi{y,x,ip). 

Therefore, we have 

Ho{-{Fo + eFi),y,x) + eHi{-{Fo + eFi),ip,y,x) = ho. 
Then let £ = : Ho{-Fo, y, x) = ho- 

So 

Ho{-{Fo + eFi),y, x) + eH^{-{Fo + eFi), tp, y, x) = Ho{-Fo, y, x), 
Ho{-{Fo + eFi),y,x) - Ho{-Fo,y,x) = -eHi{-{Fo + eFi),ip,y,x). 

Applying the intermediate value theorem here, we get that there exist some value /, such 
that 



dHo{I,y,x) 
dl 

so Fi{y,x,(p) = 
We know from assumptions that 



-sFi) = -6Hi{-{Fo + eF,),ip,y,x) 
Hi{I,(f,y,x) 



dHo{I,y,x) 



dl 



dl 



> const 



and \H,{I,^,Y{t,Io),X{t,Io))\< 



JuJo{Io,Y{n,Io),X{n,Io))dn 



2+u 



Therefore, 



|Fi(F(e,/o),X(e,/o),9^)| < const-|i/i(/,9?,F(r,/o),X(r,/o))| 

const 



< 



1 + 



2+u ■ 



□ 
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Corollary 5.1. For if e = - 6^, ^ e D/: = - 6^, 



aFi(F(e,/o),X(e,/o),<^,/io) . const 




1 + |^|2+.. 



where 6^ is a positive constant. 

Now let us draw the phase portrait of Fo{y,x) = —I: 



Figure 13 



and introduce new coordinates rj). Here rj is exactly the value of Fq on level lines Fq = r], 
and ^ is the slow time of motion. The initial section is X(0, J), Y{0,I) parametrized by /. 
Then construct a transformation (^,77) ^ {x,y). 

Lemma 5.3. In the domain x D^i, where is what mentioned above and is a 
neighbourhood of the real point —I^, the mapping r/) i-> {x,y) is analytic and canonical. 

Proof. With map {C,, rj) h-> (x, y), consider x and y as functions of ^, rj: 

x = x{^,r]), y = y{^,r]). 

From the motion above, we have obtained the relations: 

dx OFq dy dFo 
dy ' (9^ dx ' 

Let us differentiate both sides of the Hamiltonian Fq{ii,x, Hq) = rj with respect to rj: 

dFo dy ^ OFq dx ^ 
dy drj dx drj 
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Replacing with corresponding derivatives, we get 

dx dy dy dx 
dt] drj 

That is indeed the Jacobi determinant: 



1. 



D{x,y) 



det 



Therefore, (x, y) (^, 77) is canonical transformation. 



dx dx 
drj 



dy dy 
d^ drj 



1. 



□ 



After this transformation, the new Hamiltonian has the form 

-Hv, ^,^) =V + £Gi{r], ^, (p). 
Function Gi is analytic and has the following estimate: 



1^1(77,^,^)1 < 



const 



Also from Cauchy estimate [5], for r/ G -D,? = -D,, — 5,,, E = — 6^, ip G D^, where Sr, 
and 6^ are positive constants, we obtain 



d(p 

dGi{r],^,p) 



drj 

dGi{r],^,p) 



< 



< 



< 



const 

const 
const 



Because of the above estimate, we can take another similar isoenergetic reduction through 
expressing rj, via /, p, ^ and constant Hq: 

Here {—rj) is the new Hamiltonian, and ^ is the new slow time. / and p are conjugate 
variables. 
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From implicit function theorem, we can obtain the similar conclusions that function Ki is 
analytic in Dj x x and estimates of Ki and its derivatives have the following forms: 

I / r const 

\Ki{I,ip,0\ < 



1 + 



2+u ' 



dl 



dip 

dK,{I,ip,0 



< 



< 



< 



const 
const 
const 



in [Dj — 6i) X D^p x D^. Here Dj is the same as the domain — -D,; 



In this latest Hamiltonian, Kq = I and the frequency is 



dl 



1. The level lines 



Im J d^i = Im ^ = S = const, < |5| < 7 


also lie in the domain and have a positive distance from the boundary of D^. The constant 
7 is the same as introduced in Assumption 3°. 

Denoting the new time as = e~^^, we can write differential equations of motion: 

dl __ dKi _ dKi 

Now consider the exact solution !{'&), ^{d) with real initial conditions /(O), </9(0) at "i^o = 
e~^^o = 0. We obtain the relation which the solution satisfies as follows: 

J 



From the Cauchy criterion [H], we can easily prove the existence of limiting values of /(■i?) as 
— > ±00. Moreover, we also know differential equations of motion: 

dip _ ^ dKi 
dl? dip 
, dip dHo , dH, 
d^ = ^+^^- 

It is evident that ip — ?■ ±00 as — ±00. Meanwhile, t — )■ ±00. Therefore, the limits when 
t — 7- ±00 are equal to those when — )■ ±00. Let us define values 



/± = lim I{t), A/ = /+ - /_ 

t— >-±oo 
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Now the system under consideration is described by the Hamiltonian 



where / and ip are conjugate canonical variables. The action-angle variables /, (p of the 
unperturbed problem are defined in region Dj x D^. C, is the slow time variable defined in 
region and time = e~^C,- The function Ki is an analytic, 27r-periodic in ip function. 
The frequency of unperturbed motion is 1. We have proved the following properties: 

1°. The functions Kq = /, Ki can be analytically extended into a complex domain Dj x 
X D^, where Dj is some neighbourhood of a given real point is a strip of some 

fixed width about the real axis, and is a strip 



|Im^| < (Ti + 5i|Re^|, |Re^|>C5r 
|Im < (72, |Re el < CsF 



where cxi + 5i • c^T = a2- The function Ki satisfies the estimate 

Here a, 5, C, v are positive constant. 
2°. The level lines 



Imy d^ = Im ^ = S = const, < < 7 


lie in the domain and have a positive distance from the boundary of D^. 

Consider a solution I{'d), ip{'&) of the Hamiltonian system with real initial conditions /(O), 
y,(0) at^ = ^o = e-^^o = 0. For 

J± = lim A/ = /+-/_, 

i9— >-±oo 

referring to the paper [2], we have the estimate of AI that is exponentially small: 

A/ = 0(e-?). 



Actually, all variables /, ip, y and x mentioned before should contain the bar symbol " 
We obtain in fact the conclusion that 

A/ = /+-/_ = 0{e~i), 7 = const > 0. 



Now return the bar. Recall that from the original Hamiltonian E{p, q, y, x), we transform to 
two forms of Hamiltonian. With generating function S{I, q, y, x) where y, x are considered 
as parameters, we obtain Ho{I,y,x) and with e~^yx + S{I,q,y,x) we obtain H(){I,y,x) + 
eHi{I,(p,y,x). Since the limiting values of I exist, now let us prove the lemma that the 
limiting value of I is equal to that of / and thus exists. 
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Lemma 5.4. 



lim I{t) = lim I{t). 

4— >itoo t— >itoo 



Proof. With generating function S{p,q,y,x), the Hamiltonian E{p,q,y,x) is canonically 
transformed to Hq(I, y, x) and p = - ' ' ' 



dq 



As we know 



dS 
dx 



dE 
dx 



dx J ^' 



so 



dS dE 

— > 0, since )■ 0, as x ^ ±oo. 

ox ax 



On the other hand, with generating function e ^yx+S{I, q, y, x), the Hamiltonian E{p, q, y, x) 
is transformed to Ho{I, y, x) + eHi{I, (p, y, x) and p = ^' '^'^ . Also we know 



dq 



dHn 



dx 



dS dE 
5 dy dy 



dS 
~dx' 



dH dS 

where x and y are some intermediate values. As — )■ 0, — — )■ 0, thus -f/i — t- 0, as 

dx dx 

X —J- ±oo. 

fdHo\^ /dHo\'2 dHo _ , , 

From assumption + > const, and -—z )■ 0, as a; — )■ ±oo, we know that 

\ dy / \ dx J dx 



dHo 



> const for big Thus, when > cq, 
dy 



dy 

constants. From y 



dHo 



dy 



> Cj^, where cq, cy are positive 



'dHo , dH, 



X 



dx (dHo , dH,\ ^, , 

e I ^17^ + £^17^ I , the phase 



dt \ dx dx j' dt \ dy dy 

portrait of Hq + eHi = const in {x, y) plane will close to that of the adiabatic approximation 
Ho = const for |a;| < Cq + 1. 



-C6 



Figure 14 



Thus for some time moments t+ and t_, we have 
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Then for t > t+, we have 



X = e 



Therefore, x{t) — >■ +00 as t ^ +00. Similarly, x{t) — )• —00 as t — > —00. 
From 

_ dS{I,q,y,x) _ _ _ . 

y = y + e = y + £V{i, ^, y, x), 

taking limiting values on both sides as i — > ±00, we obtain y — y — > 0. 
We also know that 

P = P{I, V, y, x) = P{I, (fi, y, x), 
q = Q{I, (p, y, x) = Q{I, (fi, y, x). 

As y — y — )■ 0, it is evident that / — / — >-0, </? — >-Oasi— >■ ±00. 
Therefore, 

lim I{t) = lim I{t). 

t—^±oo t— >-±oo 

□ 

Thus / and / have the same limit, and A7 = A/. Therefore, we have the following corollary: 
Corollary 5.2. The estimate 

A7 = 0(e-?) 

is valid. 
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